function hanalysis(filename)

  load(filename, 'params', 'tOut', 'uOut');
  params
  params.plotHydro = true;

  % initialize derivative kernels
  hinit(params);

  Nx = params.Nx;
  Ny = params.Ny;
  Nt = Nx.*Ny

  time_steps = 1:params.Tsteps;
  for t=time_steps
    t = tOut(t);
    u = squeeze(uOut(t,:));


    T  = reshape(u(0*Nt+1 : 1*Nt), Nx, Ny);
    ux = reshape(u(1*Nt+1 : 2*Nt), Nx, Ny);
    uy = reshape(u(2*Nt+1 : 3*Nt), Nx, Ny);
    ut = reshape(u(3*Nt+1 : 4*Nt), Nx, Ny);

    Tdx  = dx(T);
    Tdy  = dy(T);
    uxdx = dx(ux);
    uxdy = dy(ux);
    uydx = dx(uy);
    uydy = dy(uy);



    hplot(reshape(tOut(t), 1,1), reshape(uOut(t,:), Nt*4, 1),'', params);
    pause(0.01);
  end
end


















